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Abstract 

We  address  the  multicast  problem  in  cognitive  radio  networks,  where  secondary  users  exploit 
channels  temporarily  unused  by  primary  users  (i.e.,  spectrum  opportunities).  The  existence  of  a  commu¬ 
nication  link  between  two  secondary  users  depends  not  only  on  the  transmission  power  of  the  secondary 
transmitter  and  the  distance  between  these  two  users,  but  also  on  the  occurrence  of  spectrum  opportu¬ 
nities.  This  dependency  on  the  occurrence  of  spectrum  opportunities  complicates  the  construction  of  an 
efficient  multicast  tree  in  cognitive  radio  networks.  By  taking  into  account  this  dependency,  we  propose 
a  low-complexity  approximation  algorithm  with  bounded  performance  guarantee  for  constructing  the 
minimum-energy  multicast  tree,  which  transforms  the  multicast  problem  into  a  directed  Steiner  tree 
problem.  We  also  demonstrate  this  dependency  by  studying  the  impact  of  the  traffic  load  of  the  primary 
network  on  the  minimum-energy  multicast  tree. 


I.  Introduction 

Multicast  can  provide  better  support  for  one  to  many  communications  than  unicast  or  broad¬ 
cast  [1],  and  it  thus  has  many  potential  applications  in  both  civil  and  military  domains,  e.g. 
streaming  media,  Internet  television,  and  delivery  of  situational  awareness  information  and  com¬ 
mands  on  the  battlefield.  One  of  the  most  significant  problems  in  implementing  the  multicast  is  to 
construct  an  energy-efficient  multicast  tree,  where  several  heuristic  or  approximation  algorithms 
have  been  proposed  in  [2-4]  for  the  conventional  wireless  networks. 

Multicast  in  cognitive  radio  (CR)  networks  has  received  little  attention.  In  a  CR  network, 
secondary  users  identify  and  exploit  channels  temporarily  and  locally  unused  by  primary  users 
without  causing  unacceptable  interference  to  primary  users  [5].  When  we  calculate  the  energy 
consumption  for  the  multicast  trees  in  CR  networks,  we  thus  need  to  consider  the  energy  used 
for  sensing  the  availability  of  the  channel  as  well  as  the  energy  used  for  transmissions.  The 
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proportion  of  the  sensing  energy  is  significant  especially  when  spectrum  opportunities  occur 
infrequently.  In  this  case,  much  more  energy  is  consumed  for  sensing  the  channel  before  each 
successful  transmission. 

Another  elusive  twist  in  constructing  an  efficient  multicast  tree  for  CR  networks  is  that  the 
occurrence  of  spectrum  opportunities  is  affected  by  the  transmission  power  of  secondary  users  [6]. 
If  a  secondary  user  uses  a  high  transmission  power  to  reach  a  relatively  large  number  of  multicast 
nodes,  it  must  wait  for  the  opportunity  that  no  primary  receiver  is  active  within  its  relatively 
large  interference  region,  which  happens  less  often.  As  discussed  above,  the  secondary  user  will 
spend  more  sensing  energy  for  the  emergence  of  one  opportunity.  If,  on  the  other  hand,  it  uses 
a  low  power,  it  can  only  reach  a  small  number  of  multicast  nodes,  and  more  transmissions,  each 
relying  on  its  own  opportunities  to  emerge,  are  needed  to  reach  all  the  multicast  nodes.  This 
tradeoff  in  choosing  the  secondary  users’  transmission  power  further  complicates  the  construction 
of  the  minimum-energy  multicast  tree  in  a  CR  network. 

To  summarize,  the  construction  of  the  minimum-energy  multicast  tree  in  a  CR  network  depends 
not  only  on  the  topology  of  the  secondary  network  which  is  essentially  determined  by  the 
transmission  powers  of  the  secondary  users,  but  also  on  the  occurrence  of  spectrum  opportunities 
which  is  determined  by  the  transmission  powers  of  the  secondary  users  and  the  traffic  load  of 
the  primary  network  [6].  As  a  consequence,  the  minimum-energy  multicast  tree  under  a  low 
primary  traffic  load  may  not  be  optimal  under  a  high  primary  traffic  load. 

By  considering  the  impact  of  the  occurrence  of  spectrum  opportunities,  we  propose  an  approx¬ 
imation  algorithm  for  constructing  the  minimum-energy  multicast  tree  of  a  CR  network  in  this 
paper.  The  basic  idea  of  this  algorithm  is  to  formulate  the  multicast  problem  as  a  directed  Steiner 
tree  problem  and  then  apply  an  approximation  algorithm  for  the  directed  Steiner  tree  problem. 
The  algorithm  delivers  a  solution  with  a  bounded  approximation  ratio  0(log2  \D\),  where  \D\  is 
the  number  of  destinations.  Let  n  be  the  number  of  secondary  users  and  k  their  maximal  degree1, 
then  the  running  time  of  the  approximation  algorithm  is  given  by  O  ((2/,'n)log  /  ;  _D|21og  D  +  kn!2k  j 
We  also  demonstrate  the  dependency  of  the  minimum-energy  multicast  tree  on  the  occurrence 
of  spectrum  opportunities  by  studying  the  impact  of  the  traffic  load  of  the  primary  network. 

1  The  degree  of  a  secondary  user  is  defined  as  the  number  of  secondary  users  which  are  within  its  maximal  transmission  range. 
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II.  Network  Model 

We  consider  a  secondary  network  overlaid  with  a  Poisson  distributed  primary  network,  where 
both  networks  adopt  a  slotted  transmission  structure.  Assume  that  the  primary  users  are  mobile 
and  the  relative  positions  of  the  secondary  users  are  static. 

A.  Primary  Network 

At  the  beginning  of  each  slot,  the  primary  transmitters  are  distributed  according  to  a  two- 
dimensional  Poisson  point  process  XPT  with  density  XPT.  To  each  primary  transmitter,  its 
receiver  is  uniformly  distributed  within  its  transmission  range  Rp.  Here  we  have  assumed  that  all 
the  primary  transmitters  use  the  same  transmission  power  and  the  transmitted  signals  undergo  an 
isotropic  path  loss.  Based  on  the  displacement  theorem  [7,  Chapter  5],  it  is  easy  to  see  that  the 
primary  receivers  form  another  two-dimensional  Poisson  point  process  XPR  with  density  XPT. 
Note  that  the  two  Poisson  processes  XPT  and  XPR  are  correlated. 

B.  Secondary  Network 

In  contrast  to  the  case  in  a  conventional  ad  hoc  network,  whether  the  communication  is 
successful  between  two  secondary  users  depends  not  only  on  the  transmission  power  of  the 
secondary  transmitter  and  the  distance  between  these  two  users,  but  also  on  the  availability  of 
the  communication  channel  ( i.e .,  the  presence  of  a  spectrum  opportunity).  The  latter  is  determined 
by  the  transmitting  and  receiving  activities  in  the  primary  network  as  discussed  below. 

The  notion  of  spectrum  opportunity  in  a  multicast  setting  is  open  for  interpretation.  Its  basic 
composition,  however,  roots  in  the  definition  of  spectrum  opportunity  in  unicast  [8].  Let  A  be 
the  secondary  transmitter  and  B  its  receiver  (see  Fig.  1).  Under  the  disk  signal  propagation  and 
interference  model,  there  exists  an  opportunity  from  A  and  B  if  the  transmission  from  A  does 
not  interfere  with  nearby  primary  receivers  in  the  solid  circle,  and  the  reception  at  B  is  not 
affected  by  nearby  primary  transmitters  in  the  dashed  circle.  The  radius  ry  of  the  solid  circle 
at  A,  referred  to  as  the  interference  range  of  the  secondary  users,  depends  on  the  transmission 
power  of  A  and  the  interference  tolerance  of  the  primary  receivers,  whereas  the  radius  Rj  of  the 
dashed  circle  (the  interference  range  of  the  primary  users)  depends  on  the  transmission  power 
of  the  primary  users  and  the  interference  tolerance  of  B. 
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Fig.  I.  Definition  of  spectrum  opportunity. 


It  is  clear  from  the  above  discussion  that  spectrum  opportunities  depend  on  both  transmitting 
and  receiving  activities  of  the  primary  users.  Furthermore,  spectrum  opportunities  are  asymmetric. 
Specifically,  a  channel  that  is  an  opportunity  when  A  is  the  transmitter  and  B  the  receiver  may 
not  be  an  opportunity  when  B  is  the  transmitter  and  A  the  receiver. 


III.  Minimum-Energy  Multicast  Tree 


Let  /i  be  the  multicast  source  and  D  the  set  of  destinations.  Let  T  be  a  multicast  tree  rooted 
at  /i  in  the  secondary  network,  Vr  the  node  set  of  T,  and  Vtl  ^  V  the  set  of  leaf  nodes,  then 
the  objective  and  the  constraint  of  the  minimum-energy  multicast  tree  problem  is  given  by 


arg  mm 

T 


s.t.  etx(y )  <  ejuty)  for  all  u  e  Vr,  /i  €  Vr,  and  D  C  Vr, 


E  ^  e{u) 
jsGVt—Vtl 

where  e(u)  denotes  the  energy  used  by  node  v  consisting  of  the  sensing  energy  and  the  trans¬ 
mission  energy,  etx{y)  the  energy  used  by  node  v  for  each  transmission,  and  cm(v)  the  maximal 
allowable  transmission  energy  for  node  u.  The  sum  of  e(u)  does  not  include  Vtl  because  the 
leaf  nodes  in  the  multicast  tree  neither  transmit  nor  detect  the  opportunity. 

In  general,  the  solution  to  the  minimum-energy  multicast  tree  should  specify  two  parameters 
for  each  secondary  node2:  the  transmission  power  and  the  group  of  its  intended  receivers.  But 


:If  a  secondary  node  is  not  included  in  the  multicast  tree,  then  its  transmission  power  is  zero  and  the  group  of  its  intended 
receivers  is  an  empty  set. 
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it  can  be  shown  in  the  proof  of  Theorem  1  that  given  the  group  of  the  intended  receivers,  the 
optimal  transmission  power  is  determined  by  the  distance  to  the  farthest  intended  receiver.  It 
follows  that  the  minimum-energy  multicast  tree  is  fully  determined  by  the  group  of  the  intended 
receivers  for  each  secondary  node. 

The  approximation  algorithm  for  the  minimum-energy  multicast  tree  is  based  on  an  auxiliary 
directed  graph  with  edge  weights.  We  will  address  the  construction  of  the  auxiliary  graph  G  and 
the  calculation  of  the  edge  weights  in  the  following  two  subsections. 

A.  Construction  of  the  Auxiliary  Graph  G 

The  auxiliary  graph  G  is  constructed  in  the  following  way.  For  each  secondary  node  p  with 
degree  k,  we  add  2k  —  1  multicast  nodes  u,  (1  <  i  <  2k  —  1)  which  represent  all  possible  groups 
of  its  intended  receivers.  Then  we  connect  the  secondary  node  p  with  each  of  its  2k  —  1  multicast 
nodes  by  a  directed  edge  and  assign  a  positive  weight  to  each  edge  (see  Fig.  2).  For  the  directed 
edge  (p,  z/j),  the  assigned  weight  is  equal  to  the  average  energy  required  for  transmitting 
successfully  from  p  to  the  group  of  the  intended  receivers  denoted  by  z/j.  Finally,  we  connect 
each  z/j  with  its  corresponding  intended  receivers  by  directed  edges  and  assign  zero  weight  to 
all  those  edges. 


v\ 


V2 


V 3 


Z/2fc_i 


Fig.  2.  Generation  of  multicast  nodes  for  node  \x  in  the  secondary  network. 


Based  on  the  above  construction,  we  see  that  G  is  an  edge- weighted  directed  graph,  and  it  has 
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two  types  of  nodes:  user  node  and  multicast  node.  The  user  nodes  represent  the  secondary  users, 
and  the  multicast  nodes  specify  the  transmission  power  and  the  group  of  intended  receivers  for 
the  user  nodes.  Let  Vq  be  the  vertex  set  of  G,  and  Eq  the  edge  set  of  G.  Given  the  source  //  e  V(j 
and  a  destination  set  D  C  Vc,  the  directed  Steiner  problem  A(/i.  D )  is  to  find  a  tree  T  rooted 
at  fi  which  spans  all  the  nodes  in  D  with  the  minimum  weight.  Then  we  have  the  following 
theorem  about  the  relationship  between  the  directed  Steiner  tree  in  G  and  the  minimum-energy 
multicast  tree  in  the  secondary  network. 

Theorem  1:  Given  a  multicast  session  in  the  secondary  network,  let  /i  be  the  source,  and 
D  the  destination  set.  Define  the  weight  of  a  tree  in  the  auxiliary  graph  G  as  the  total  sum 
of  the  weights  of  the  edges  in  the  tree,  and  the  weight  of  a  multicast  tree  in  the  secondary 
network  as  the  total  sum  of  the  average  energy  used  by  the  nodes  in  the  multicast  tree.  Then 
the  directed  Steiner  tree  rooted  at  //  which  spans  all  the  nodes  in  D  has  the  same  weight  as  the 
minimum-energy  multicast  tree  in  the  secondary  network. 

Proof:  Although  the  set  of  the  transmission  power  that  each  secondary  node  can  choose  is 
continuous,  it  suffices  to  consider  a  finite  subset  of  it  in  order  to  obtain  the  minimum-energy 
multicast  tree.  For  a  secondary  node  /i  in  the  minimum-energy  multicast  tree,  let  /xi,  /x 2,..., 
fik  be  the  secondary  nodes  within  the  maximum  transmission  range  of  /x  which  are  listed  in 
an  ascending  order  of  the  distance  from  /x.  Then  under  the  disk  signal  propagation  model, 
the  optimal  transmission  power  for  /x  must  belong  to  {0,  Cdf,  Cd% , ...,  Cdk},  where  C  is  a 
constant,  dj  (1  <  i  <  k)  is  the  distance  between  p  and  /i,:,  and  a  path-loss  exponent.  Let 
d0  =  0  and  dk+1  =  oo.  Suppose  that  the  optimal  transmission  power  takes  some  other  value 
p*  £  {0,  Cd",  Cd%, ...,  Cdk}  with  Cdf  <  p*  <  Cdf+1  for  some  i.  We  can  always  replace  p*  by 
Cdf  for  the  node  /x,  and  it  does  not  change  the  multicast  tree  but  reduces  the  energy  cost  of 
the  multicast  tree.  It  follows  that  the  minimum-energy  multicast  tree  is  fully  determined  by  the 
group  of  intended  receivers  for  each  secondary  node.  Based  on  the  construction  of  the  auxiliary 
graph  G,  we  thus  conclude  Theorem  1.  ■ 

Based  on  Theorem  1,  the  problem  of  the  minimum-energy  multicast  tree  is  transformed 
into  a  directed  Steiner  problem  via  the  construction  of  the  auxiliary  graph  G.  By  applying 
an  approximation  algorithm  for  the  directed  Steiner  problem  [9]  to  G  (see  Sec.  IV),  we  can 
obtain  an  approximate  solution  to  the  problem  of  the  minimum-energy  multicast  tree. 
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B.  Calculation  of  the  Edge  Weights 

Let  v  denote  the  multicast  node  which  is  connected  to  a  group  of  intended  receivers  {/xi,  /i2, /xTO} 
for  secondary  node  /i,  then  the  weight  w  of  the  directed  edge  (/x,  u)  is  the  expectation  of  the 
sum  of  sensing  energy  and  transmission  energy  until  the  last  one  in  the  group  {/xi,  /x2,  ...,/xm} 
receives  successfully  from  /x,  and  it  can  be  expressed  as 


(1) 


t= l 


where  t  is  the  slot  index,  T  a  stopping  time  which  is  the  number  of  slots  for  /x  transmitting 
successfully  to  all  the  receivers  that  are  connected  to  v,  es  the  sensing  energy  per  slot,  etx  the 
energy  used  for  each  transmission,  and  I tx{t)  an  indicator  function  which  is  1  if  /x  transmits  in 
slot  t  and  0  otherwise. 

Consider  the  following  two  transmission  schemes  for  the  secondary  node  /x:  one  is  that  /x 
transmits  even  if  /x  only  sees  the  opportunity  from  /x  to  some  of  the  receivers  that  are  connected 
to  the  other  is  that  p  does  not  transmit  until  /x  sees  the  opportunity  from  /x  to  all  the  receivers 
that  are  connected  to  v.  If  /x  uses  the  first  transmission  scheme  (referred  to  as  sequential  scheme), 
it  may  possibly  transmit  more  than  once,  where  each  time  it  transmits  to  disjoint  subgroups  of 
receivers  that  are  connected  to  u.  Although  the  calculation  of  (1)  can  be  formulated  into  solving 
the  absorbing  time  of  a  Markov  chain  in  this  case,  the  obtained  expression  is  complicated  and 
difficult  to  evaluate. 

If  /x  adopts  the  second  transmission  scheme  (referred  to  as  simultaneous  scheme),  /x  transmits 
only  once  to  all  the  receivers  that  are  connected  to  v  when  /x  sees  the  opportunity  from  /x  to  all 
the  receivers.  Then  in  this  case,  the  edge  weight  w  can  be  rewritten  as 


w  =  E  [T)es  +  etx, 


where  T  is  the  first  slot  that  /x  sees  the  opportunity  from  /x  to  all  the  receivers  that  are  connected 
to  u.  Due  to  the  i.i.d.  distribution  of  the  primary  network  over  slots,  T  is  obviously  a  geometric 
random  variable  with  parameter  p0,  where  p0  is  the  probability  of  having  an  opportunity  from  /x 
to  all  the  receivers  {/xi,  /x2, ...,  /xm}  at  any  given  time.  Under  the  disk  signal  propagation  model, 
we  thus  have  that 


w 


t's  / Po  T  &tx 
es/p0  +  Cda, 
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where  C  is  a  constant,  d  the  distance  between  p  and  its  farthest  intended  receiver,  and  a  the 
path-loss  exponent. 

The  following  proposition  gives  the  expression  for  p0. 

Proposition  1:  Let  A pt  be  the  density  of  the  primary  transmitters.  Let  Rp  and  Rj  denote 
the  transmission  range  and  the  interference  range,  respectively,  of  the  primary  users.  Given 
a  secondary  user  p  and  a  group  of  its  intended  receivers  { / /  , ,  p2,  ••••  /Jm],  let  ry  denote  the 
interference  range3  of  /1  when  its  transmission  reaches  all  the  intended  receivers.  Choosing  p  as 
the  origin  of  the  polar  coordinate  system,  we  have  that  the  probability  p0  of  having  an  opportunity 
from  p  to  all  its  intended  receivers  is  given  by 


Si(r,  Rp,rf)  is  the  common  area  of  two  circles  with  radii  Rp  and  ry  and  centered  r  apart,  and 
S(j  is  the  union  of  the  circles  with  radii  Rj  centered  at  each  intended  receiver  pt  (1  <  1  <  rri). 

Furthermore,  we  have 

exp  l-XpT^r]  +  Su)]  <  p0  <  min  {exp  (-XPT7t rj)  ,  exp  (-XpTSu)}  ,  (3) 

where  Su  is  the  area  of  the  region  Su. 

Proof:  From  the  definition  of  unicast  spectrum  opportunity  given  in  Sec.  II-B,  we  know 
that  a  spectrum  opportunity  from  p  to  all  its  intended  receivers  {pi,  p2: ....  pm}  occurs  if  and 
only  if  there  are  no  primary  receivers  within  distance  ry  of  p  and  no  primary  transmitters  within 
distance  Rj  of  each  receiver  pt  (1  <i  <m). 

Let  I (p,  d,  rx/tx)  denote  the  event  that  there  exists  primary  receivers/transmitters  within  dis¬ 
tance  d  of  a  secondary  user  p.  Let  ![(//,  d,  rx/tx)  denote  the  complement  of  I(/r,  d,  rx/tx).  Then 


3Since  the  minimum  transmission  power  for  successful  reception  is,  in  general,  higher  than  the  maximum  allowable  interference 
power,  it  follows  that  the  transmission  range  rp  of  p  is  smaller  than  its  interference  range  r/,  where  rp  is  set  to  be  the  farthest 
distance  between  /i  and  its  intended  receivers.  Furthermore,  under  the  disk  signal  propagation  and  interference  model,  we  have 

rp  =  fin  (0  <  /3  <  1). 
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the  probability  p0  of  having  an  opportunity  from  /i  to  {/ii,  p2,  ■  ■■,  pm}  is  given  by 


Po  =  Pr{%,r7,rx)  nl(/ii,f?7,tx)  n  ■  ■  •  n%m,  f?/,tx)} 


=  Pr{I(/ii7JR/,tx)n---nl(/im,  f?/,tx)|  I(/i,r/,rx)}Pr{I(/i,r/,rx)}.  (4) 

Since  the  primary  receivers  admit  a  Poisson  point  process  with  density  A pt,  it  follows  that 

Pr{%,ry,rx)}  =  exp  {-\PTTir])  .  (5) 


Given  I(yU,  77,  rx),  i.e.,  there  are  no  primary  receivers  within  distance  ry  of  p,  based  on 
Independent  Thinning  Theorem  [7,  Chapter  5],  we  have  that  the  primary  transmitters  form  an 
inhomogeneous  Poisson  point  process  with  density  specified  by 

Sj(r,  Rp,  r j)' 


ApT(r)  —  A  pt 


1  - 


nR2p 


Thus, 


Pr{II(/ii5  Rh  tx)  n  ■  ■  •  nl(/im,i?/,tx)|  I(/i,r7,rx)}  =  exp  |  -  J  J  A/PT(r)rdrd6'  |  .  (6) 

Su 

Plugging  (5,  6)  into  (4)  yields  (2). 

From  (4),  we  can  easily  see  that 


Po  <  min{Pr{I(/i,r7,rx)},  Pr{I(^1,i?/,tx)n---nl(^m,  i?j,tx)}} 
=  min{exp  (-A ,  exp  (-\PTSu)}. 

On  the  other  hand,  it  follows  from  (6)  that 


Pr{I(/il7  Rp  tx)  fl  •  ■  ■  n  Rp  tx) I  I(n,  rp  rx)}  >  Pr{I(^l7  RI?  tx)  (T  ■  ■  ■  (T  I(/xm,  Ri,  tx)}. 


By  considering  (4),  we  have  that 

Po  >  Pr{I(/r,  vy,  rx)}Pr{I(/ii,  Rj,  tx)  fl  ■  ■  •  fl  Ri,  tx)} 

=  exp  [-XpT^r]  +  Sv)]  . 

Thus,  the  inequality  (3)  is  shown.  ■ 

From  (2),  we  can  see  that  the  computation  of  p0  requires  a  double  integral  over  an  irregular 
region  which  is  the  union  of  several  circles.  Due  to  the  complexity  of  the  computation,  we  use 
the  lower  bound  on  p0  given  in  (3)  as  an  approximation  of  p0  to  obtain  the  edge  weight  w,  i.e., 

6 s 

W  ~ 


exp  [-APT(7rrf  +  S}/)] 


+  Cda. 


(7) 
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IV.  Approximation  Algorithm  for  Directed  Steiner  Trees 


In  Sec.  Ill,  we  transform  the  problem  of  the  minimum-energy  multicast  tree  into  the  directed 
Steiner  problem  by  deriving  a  directed  graph  G  with  weighted  edges  from  the  secondary  network. 
In  this  section,  we  present  an  approximation  algorithm  for  the  directed  Steiner  problem  [9]. 

Consider  a  directed  graph  G  =  (VG,  EG)  where  VG  is  the  vertex  set,  EG  is  the  edge  set,  and 
each  edge  e  e  EG  is  associated  with  a  nonnegative  weight  w(e).  Given  a  source  /i  <E  VG  and 
a  destination  set  D  C  VG,  we  recall  that  the  directed  Steiner  problem  A(ji.  D )  is  to  find  a  tree 
T  rooted  at  fi  which  spans  all  the  nodes  in  D  with  the  minimum  weight.  Let  k  —  \D\  be  the 
number  of  destinations.  A  slightly  more  general  version  of  the  directed  Steiner  problem  is  to 
find  a  tree  T  rooted  at  //  which  spans  at  least  m  nodes  in  D  with  the  minimum  weight  for  some 
m  <  k,  which  is  denoted  by  A(n,  m,  D). 

1)  Preliminary  Results:  Assume  that  the  directed  graph  G  is  strongly  connected,  i.e.,  for 
every  pair  of  nodes  fi  and  //,  there  exist  at  least  one  directed  path  in  G  from  //  to  u  and  at 
least  one  directed  path  in  G  from  v  to  ji.  Then  based  on  G,  we  construct  a  complete  graph  G 
which  has  the  same  vertex  set  VG  as  G,  and  we  let  the  weight  v)  of  each  edge  (/x,  v)  in 
G  equal  the  shortest  path  distance  from  //  to  u  in  G.  Although  G  contains  more  edges  than  G, 
the  solution  T  to  the  directed  Steiner  problem  A(fi,  D)  directly  leads  to  the  solution  T  to  the 
directed  Steiner  problem  *4.(/x,  D )  if  we  replace  the  edges  of  T  that  do  not  exist  in  G  by  the 
corresponding  shortest  paths  in  G. 

Let  if;  (T)  denote  the  weight  of  a  tree  T  which  is  the  total  sum  of  the  weights  of  the  edges 
in  T,  and  k(T)  the  number  of  destinations  contained  in  T,  i.e.,  k{T )  =  j Vr  fl  D |.  Then  we 
define  the  density  d(T)  of  tree  T  as  the  ratio  of  the  weight  of  T  to  the  number  of  destinations 
contained  in  T,  i.e., 


M  T)k<P 
1  k(T) 

A  tree  is  said  to  be  an  /-level  tree  if  no  leaf  is  more  than  l  edges  away  from  the  root.  It  is 
shown  in  [10]  that  for  all  Z  >  1  there  exists  an  /-level  tree  that  provides  a  A-l///  approximation 
to  the  general  directed  Steiner  problem  A(ii,m,D),  where  m  <  k  is  the  minimum  number  of 
destinations  required  to  be  contained  in  the  /-level  tree.  A  trivial  approximation  when  /  =  1  is 
the  combination  of  the  edges  of  G  from  the  root  //  to  the  m  destinations  which  are  closer  to  // 
than  the  other  k  —  m  destinations  (in  terms  of  the  shortest  path  distance).  This  approximation  to 
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A(n,  m,  D )  gives  an  approximation  to  A(n,  m,  D )  if  we  change  the  edges  of  G  that  do  not  exist 
in  G  into  the  corresponding  shortest  paths  in  G .  Let  Tx  denote  the  above  approximate  solution 
to  A(fx,  rn.  D )  and  T*  the  directed  Steiner  tree  for  A(ji,  m,  D )  (i.e.,  the  optimal  solution).  It  can 
be  easily  shown  that  the 


uATO 

w(T*) 


<  m. 


2 )  Algorithm:  This  approximation  algorithm  is  a  special  case  of  the  approximation  algorithm 
for  the  general  directed  Steiner  tree  problem  A(n,m,  D)  [9].  The  basic  idea  is  to  recursively 
obtain  a  [log fc] -level  tree  T  which  provides  a  log/c(logfc  —  l)/c1//logfc  approximation  to  the 
directed  Steiner  problem  A{n,  D )  and  then  transform  T  into  a  tree  T  in  the  original  directed 
graph  G. 

Approximate  Algorithm  for  Directed  Steiner  Problem  A(n,  D) 


INPUT: 


(i)  An  edge- weighted  directed  graph  G  =  (VG,EG),  where  each  edge  e  6  EG  is  associated 
with  a  nonnegative  weight  w(e); 

(ii)  A  source  /i  <E  VG  and  a  destination  set  D  C  VG  with  size  k  =  \D\. 

OUTPUT:  A  tree  T  rooted  at  /i  that  spans  all  the  nodes  in  I). 

STEP  0  If  not  all  the  destinations  in  D  are  reachable  from  /i,  then  RETURN  0. 

STEP  1  Derive  a  complete  graph  G  from  G.  Specifically,  for  every  pair  of  nodes  //  and  v  in  G,  if 
there  exists  a  directed  edge  (fx,  v)  from  n  to  v  in  G,  i.e.,  (/i,  v)  e  E,  we  keep  the  edge 
(/i,  v)  along  with  its  weight  in  G\  if  not,  i.e.,  (/i,  v)  0  E,  we  add  an  edge  (/i,  v)  in  G  and 
assign  the  shortest  path  distance  from  fi  to  v  in  G  to  the  edge  as  its  weight  in  G. 

STEP  2  Invoke  the  function  k,  D )  which  returns  an  I- level  tree  T  as  an  approximate  solution 
to  the  directed  Steiner  problem  A(fx,D),  where 

I  =  [log  A:]. 

Given  m  <  \D\,  the  function  F,(/l  rn,  D )  for  i  >  1  is  defined  recursively  as  follows: 
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k  <—  \D\  (the  number  of  nodes  in  D). 

IF  i  =  1 

T  <—  the  combination  of  the  edges  of  G  from  /x  to  the  m  nodes  in  D 
which  are  closer  to  /i  than  the  other  k  —  m  nodes  in  D. 

RETURN  f. 

ELSE 

T<-0. 

Die  ft  D  (the  set  of  destinations  that  are  not  contained  in  T). 
mieft  <—  m  (the  number  of  destinations  that  are  needed  by  T). 

WHILE  mieft  >  0 

Dbest  *  0  and  d(Tbest)  oo. 

FOR  each  vertex  v  6  V  and  each  m'  (1  <m'<  mieft ) 

V  m',  Dieft)  U  {(/i,  u)}. 

IF  d(T')  <  d{Tbest )  (compare  the  density  of  the  two  trees) 

Tbest  <—  T'. 

T<-TU  Tbest  (combine  T  and  Tbest  into  one  tree). 

TTlleft  TYlleft  | Dieft  H  V (Tbest)  |. 

Dieft  Dieft  1  ( Dbest )  (U  {Tbest)  is  the  vertex  set  of  Tbest)’ 

RETURN  T. 

Notice  that  all  the  operations  of  the  function  Fj(/x,  m,  D)  are  performed  on  the  complete 
graph  G,  not  the  original  directed  graph  G. 

STEP  3  Transform  the  tree  T  in  G  into  a  tree  T  in  G  by  changing  the  edges  of  T  that  do  not  exist 
in  G  into  the  corresponding  shortest  paths  in  G.  The  detailed  procedure  is  shown  as  below. 
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T<-0. 

FOR  each  edge  e  <E  E(T)  ( E(T )  is  the  edge  set  of  tree  T ) 
IF  e  e  f?(G)  ( E(G )  is  the  edge  set  of  graph  G) 

T  <-  T  U  e. 

ELSE 

Le  <—  the  shortest  path  in  G  corresponding  to  e. 
T^TULe. 

RETURN  T. 

This  is  the  end  of  this  algorithm. 


Remark.  If  T  is  an  optimal  solution  to  the  directed  Steiner  problem  A(fi,  D ),  then  after  we 
replace  the  edges  of  T  that  do  not  exist  in  G  by  the  corresponding  shortest  paths  in  G,  we 
obtain  an  optimal  solution  T  to  the  directed  Steiner  problem  A(/jl,  D).  It  is  easy  to  see  that  T 
does  not  have  duplicate  edges  of  G,  otherwise  removing  these  additional  edges  yields  a  tree  T' 
with  a  smaller  weight  than  T  which  leads  to  a  contradiction.  But  since  T  is  only  an  approximate 
solution  to  the  directed  Steiner  problem  A(n,  D ),  it  is  possible  that  during  step  3,  before  adding 
edge  e  or  the  shortest  path  Le  corresponding  to  edge  e  to  the  tree  T,  either  the  edge  e  or  some 
edge  on  the  shortest  path  Le  is  already  included  in  the  tree  T.  In  this  case,  this  duplicate  edge 
is  not  added  to  the  tree  T  after  the  combination  of  T  and  e  or  Le.  It  implies  that 

w(T)  <  w(T). 


On  the  other  hand,  it  follows  from  the  relation  between  the  original  directed  graph  G  and  the 
complete  graph  G  that  w(T*)  =  w(T*)  where  T*  and  T*  are  the  optimal  solutions  to  the  directed 
Steiner  problems  A(fi,  D )  and  Aiji.  D ),  respectively.  We  thus  have  that  the  approximation  ratio 
of  T  is  bounded  above  by  the  approximation  ratio  of  T,  i.e., 


w(T)  <  w(T) 


w(T')  ~  W(T-) 

where  k  =  \D\  is  the  number  of  destinations  in  D. 


0(log2  k), 


V.  Approximation  Algorithm  for  Minimum-Energy  Multicast  Tree 

In  this  section,  we  present  the  approximation  algorithm  for  the  minimum-energy  multicast 
tree  and  analyze  its  approximation  ratio  and  time  complexity. 
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A.  Algorithm 

As  discussed  in  Sec.  Ill,  the  approximation  algorithm  for  the  minimum-energy  multicast  tree 
relies  on  the  approximation  algorithm  for  the  directed  Steiner  problem  given  in  Sec.  IV.  The 
procedure  of  this  approximation  algorithm  is  detailed  as  below. 

Approximation  Algorithm  for  Minimum-Energy  Multicast  Tree 

INPUT: 


(i)  A  secondary  network  where  V  denotes  the  set  of  nodes. 

(ii)  The  transmission  range  of  each  secondary  node  which  is  determined  by  the  maximal 
transmission  power  for  the  secondary  node. 

(iii)  A  source  fi  e  V  and  a  destination  set  D  C  V  with  size  k  —  \D\. 

OUTPUT:  A  multicast  tree  T  rooted  at  /i  that  spans  all  the  nodes  in  D. 

STEP  1  Construct  the  auxiliary  graph  G. 

1.1  Add  the  multicast  nodes  for  each  secondary  node  (see  Sec.  III-A). 

1.2  Use  (7)  to  compute  the  weight  for  each  edge  in  G.  Since  there  is  no  closed-form 
expression  for  Sjj  which  is  the  area  of  the  union  Sj-  of  several  circles,  we  resort  to 
the  Monte  Carlo  method.  Specifically,  we  generate  N  points  uniformly  distributed  in 
a  square  with  side  length  ds  which  contains  Su,  and  we  count  the  number  of  points 
that  fall  into  Su,  denoted  by  Ns.  Then  we  have 


Su 


STEP  2  Apply  the  approximation  algorithm  for  the  directed  Steiner  problem  given  in  Sec.  IV  to 
the  auxiliary  graph  G  and  obtain  an  approximate  Steiner  tree  TG  in  G. 

STEP  3  Transform  the  approximate  Steiner  tree  Tq  into  a  multicast  tree  T  in  the  secondary  network. 
The  detailed  procedure  is  shown  in  the  following  table. 
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FOR  each  node  v  e  V (TG)  (V (TG)  is  the  vertex  set  of  tree  TG) 

IF  v  e  V  ( i.e .,  v  is  a  user  node,  where  V  is  the  node  set  of  the  secondary  network) 
T  <-  T  U  v. 

ELSE  (i.e.,  v  is  a  multicast  node) 

{vx,v2,  ...,vm}  the  group  of  intended  receivers  denoted  by  v. 

T  <-  TU  ...,vm}. 

/i  <—  the  parent  of  v  in  TG. 

T  <—  T  U  (U”f1(/i,  Vi)),  where  (/x,  z/j)  is  a  directed  edge  from  /i  to  iy. 

RETURN  T. 

This  is  the  end  of  the  algorithm. 

Remark.  It  is  possible  that  for  some  user  node  v  in  the  approximate  Steiner  tree  TG,  it  has 
more  than  one  child,  that  is,  it  has  more  than  one  group  of  intended  receivers.  Take  two  groups 
G i,  G 2  of  intended  receivers  as  an  example,  and  the  case  of  more  than  two  groups  of  intended 
receivers  can  be  easily  extended.  If  G\  fl  G2  =  0,  it  implies  that  the  secondary  node  v  should 
first  multicast  to  G\  and  then  G2  (a  sequential  way),  or  multicast  to  G\  in  odd  slots  and  to  G2 
in  even  slots  (an  alternating  way).  If  Gi  fl  G2  7^  0,  then  we  can  replace  G1  by  G\  —  G 1  fl  G2 
and  do  the  multicasting  as  the  case  of  G\  fl  G2  =  0.  It  is  easy  to  see  that  this  replacement  of  G\ 
can  only  reduce  the  average  energy  cost  of  the  resulting  multicast  tree,  and  it  does  not  change 
the  resulting  multicast  tree. 

B.  Approximation  Ratio  and  Complexity 

For  the  approximation  ratio  and  the  complexity,  we  have  the  following  theorem. 

Theorem  2:  Let  T*  be  the  minimum-energy  multicast  tree  in  the  secondary  network,  and  T 
the  multicast  tree  given  by  the  approximation  algorithm  in  Sec.  V-A.  Define  the  weight  w(T) 
of  a  multicast  tree  T  as  the  total  sum  of  the  average  energy  used  by  the  nodes  in  the  multicast 
tree.  Then  we  have 

<  RoglDIK Region  -  1)101-'“*™  (8) 

where  D  is  the  set  of  destinations  and  \D\  denotes  the  size  of  D. 
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Furthermore,  the  time  complexity  of  the  approximation  algorithm  for  the  minimum-energy 
multicast  tree  is  O  ^(2fcn)1°f’ 1  |.D|21og l-°l  +  kn2k^j,  where  k  is  the  maximal  degree  of  the  sec¬ 
ondary  nodes  and  n  the  number  of  secondary  nodes. 

Proof:  Based  on  Theorem  1,  we  know  that  the  weight  w(T*)  of  the  minimum-energy 
multicast  tree  T  is  equal  to  the  weight  w(Tq)  of  the  Steiner  tree  Tf .  Moreover,  when  we 
transform  the  approximate  Steiner  tree  TG  into  the  multicast  tree  T  in  Step  3,  we  have  that 
w(T)  <  w(Tq).  Since  w(Tg)/w(Tq )  <  log  | ZT | (log  | ZT |  —  l)|i9|~loglDl  (Theorem  4  in  [9]),  it 
follows  that 

<  ^||  <  w(Tc)/w(TZ)  <  log  |D|(log  |D|  -  l)|D|-'°^. 

On  the  other  hand,  Step  1  takes  0{kn2k )  time  since  every  user  node  has  at  most  2fc  — 1  multicast 
nodes,  and  each  multicast  node  is  connected  to  at  most  k  intended  receivers.  Similarly,  Step  3 
can  be  done  in  0{kn2k)  time.  By  noticing  that  the  auxiliary  graph  G  has  0{2kn)  nodes  and 
setting  i  =  log  \D\  in  Theorem  4  [9],  we  have  that  Step  2  takes  O  ^(2 fcn)log'D'  |£)|21oglI?l  j  time. 
Thus,  the  time  complexity  of  the  approximation  algorithm  is  O  ((2l:n|1""  '  I)\iv,g  r>  +  ktV2l  \ . 


VI.  Simulation  Results 

In  this  section,  we  present  several  simulation  results.  We  place  20  secondary  nodes  uniformly 
in  a  500m  x  500m  square  such  that  they  are  connected  if  they  use  their  maximal  transmission 
power,  and  we  fix  their  positions  during  the  whole  simulation.  Other  simulation  parameters  are 
given  by:  rp  =  Rp  =  130m,  ry  =  Rj  =  144m,  etx  =  3  x  10_4d3,  and  es  =  1. 

Since  there  is  no  comparable  algorithm  for  cognitive  radio  networks,  we  choose  the  approx¬ 
imation  algorithm  proposed  in  [4]  for  conventional  wireless  networks  and  compare  its  perfor¬ 
mance  with  that  of  our  approximation  algorithm  in  cognitive  radio  networks.  The  percentage 
of  the  average  energy  saving  of  our  approximation  algorithm  under  different  traffic  load  of  the 
primary  network  (represented  by  the  density  of  the  primary  transmitters)  is  shown  in  Fig.  3. 
The  conventional  approximation  algorithm,  which  treats  the  secondary  network  the  same  as  a 
primary  network,  transforms  the  multicast  problem  into  an  undirected  node-weighted  Steiner 
tree  problem,  and  gives  a  multicast  tree  invariant  to  the  primary  traffic  load.  Conversely,  our 
approximation  algorithm  takes  into  account  the  sensing  energy,  and  produces  the  multicast  trees 
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that  adapt  to  the  primary  traffic  load.  We  can  see  that  our  approximation  algorithm  can  save  up 
to  9  ~  34%  of  the  total  energy,  and  the  energy  saving  becomes  more  substantial  as  the  primary 
traffic  load  increases. 


Fig.  3.  Percentage  of  average  energy  saving  vs.  density  of  primary  transmitters. 

We  also  study  the  impact  of  the  primary  traffic  load  on  the  minimum-energy  multicast  tree. 
Fig.  4  shows  the  average  energy  cost  of  the  sequential  and  the  simultaneous  schemes  (see 
Sec.  III-B)  for  a  fixed  multicast  tree.  We  observe  that  the  sequential  scheme  performs  worse 
when  the  primary  traffic  load  is  low,  whereas  it  performs  better  when  the  primary  traffic  load 
is  high.  Since  the  sequential  scheme  may  probably  lead  to  more  than  one  transmission  for 
each  edge  in  the  multicast  tree,  it  usually  costs  more  transmission  energy  than  the  simultaneous 
scheme.  But  on  the  other  hand,  it  does  not  need  to  wait  for  all  the  intended  receivers  seeing 
the  opportunity,  and  thus  it  costs  less  sensing  energy  than  the  simultaneous  scheme.  As  the 
primary  traffic  load  increases,  its  saving  in  the  sensing  energy  gradually  exceeds  its  extra  cost 
in  transmission  energy. 


VII.  Conclusion 

In  this  paper,  we  have  presented  a  low-complexity  approximation  algorithm  with  bounded 
performance  guarantee  for  constructing  the  minimum-energy  multicast  tree  in  cognitive  radio 
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Fig.  4.  Average  energy  cost  of  the  sequential  and  the  simultaneous  schemes  (in  logarithmic  scale)  vs.  density  of  primary 
transmitters. 


networks.  This  approximation  algorithm  takes  into  account  the  energy  used  for  sensing  the 
spectrum  opportunities,  and  its  constructed  multicast  trees  are  adaptive  to  the  traffic  load  of 
the  primary  network.  Simulation  results  have  demonstrated  its  substantial  energy  saving  energy, 
compared  with  a  well-known  conventional  approximation  algorithm.  Furthermore,  we  study  the 
impact  of  the  primary  traffic  load  on  the  minimum-energy  multicast  tree.  Simulation  results  have 
shown  that  the  simultaneous  transmission  scheme  is  more  suitable  for  light  primary  traffic  load, 
while  the  sequential  transmission  scheme  is  desired  for  heavy  primary  traffic  load. 
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